------------------------------------------------------------------------------------------------------------------------------------
      name:  plog_847
       log:  /accounts/projects/jr_ra/GRscarring/erratum/programs/analysis/fig_ur_age.log
  log type:  text
 opened on:  27 Nov 2024, 17:55:01

. *************************************************************************************************************** 
. * NG: 9/11/2019
. * fig_ur_age.do
. * Creates Unemployment and non-employment figures for young college graduates 
. * Edit: NG 01/10/20: Separate fig_ur_age in two groups and change some cosmetics aspects of the figures
. *               NG 02/7/20: Updated to a 7-month smoothing and changed to two panels 
. *.      JR 4/1/20:   Rewritten completely
. *************************************************************************************************************** 
. 
. clear

. 
. cap project, doinfo

. //cap notaproject
. if _rc==0 {
.         local pdir "`r(pdir)'"                                                      // the project's main dir.
.         local dofile "`r(dofile)'"                                                  // do-file's stub name
.         local sig {bind:{hi:[`dofile'.dta. RP : `dofile'.do, `c(current_date)']}}       // a signature in notes
.         local doasproject=1
. }

. else {
.         local pdir "~/GRscarring"
.         local dofile "fig_ur_age"
.         local doasproject=0
. }

. 
. set more off

. local rootdir "`pdir'"

. local thisdir "`pdir'"

. 
. local scratch "`pdir'/scratch"

. local rawdata "`pdir'/rawdata"

. local output "`pdir'/results"

. 
. if `doasproject'==1 {
.         project, uses("`scratch'/extractcps.dta.gz")
project GRscar_erratum > do-file uses: "/scratch/public/jr_ra/GRscarring2024/erratum/scratch/extractcps.dta.gz" filesig(1117482818:1
> 425161591)
. }

. 
. set scheme s1color

. 
. !zcat `scratch'/extractcps.dta.gz > `scratch'/extractcps.dta 


. 
. // Program to implement simple seasonal adjustment and smoothing
. cap program drop seasadj

. program define seasadj
  1.   args v
  2.   *Simple seasonal adjustment
.   areg `v' i.year, a(mon)
  3.   predict monfx, d
  4.   gen `v'_sa=`v'-monfx
  5.   drop monfx
  6. end

. 
. cap program drop smoother

. program define smoother
  1.   args v
  2.   tsset yearmo
  3.   gen     `v'_m7=100*((4*`v' + 3*L.`v' + 2*L2.`v' + 1*L3.`v'+ 3*F.`v'+ 2*F2.`v' + 1*F3.`v')/16)
  4.   replace `v'_m7=100*((4*`v' + 3*L.`v' + 2*L2.`v' + 1*L3.`v'+ 3*F.`v'+ 2*F2.`v' )/15) if F3.`v'==.
  5.   replace `v'_m7=100*((4*`v' + 3*L.`v' + 2*L2.`v' + 1*L3.`v'+ 3*F.`v' )/13) if F2.`v'==.
  6.   replace `v'_m7=100*((4*`v' + 3*L.`v' + 2*L2.`v' + 1*L3.`v')/10) if F.`v'==.
  7.   replace `v'_m7=100*((4*`v' + 3*L.`v' + 2*L2.`v' + 3*F.`v'+ 2*F2.`v' + 1*F3.`v')/15) if L3.`v'==.
  8.   replace `v'_m7=100*((4*`v' + 3*L.`v' + 3*F.`v'+ 2*F2.`v' + 1*F3.`v')/13) if L2.`v'==.
  9.   replace `v'_m7=100*((4*`v' + 3*F.`v' + 2*F2.`v' + 1*F3.`v')/10) if L.`v'==.
 10. 
.   gen     `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v' + 1*L4.`v' ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v' + 1*F4.`v')/25)
 11.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v' + 1*L4.`v' ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v'           )/24) if F4.`v'==.
 12.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v' + 1*L4.`v' ///
>                              + 4*F.`v' + 3*F2.`v'                      )/22) if F3.`v'==.
 13.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v' + 1*L4.`v' ///
>                              + 4*F.`v'                                 )/19) if F2.`v'==.
 14.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v' + 1*L4.`v' ///
>                                                                        )/15) if F.`v'==.
 15.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v' + 2*L3.`v'  ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v' + 1*F4.`v')/24) if L4.`v'==.
 16.   replace `v'_m9=100*((5*`v' + 4*L.`v' + 3*L2.`v'   ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v' + 1*F4.`v')/22) if L3.`v'==.
 17.   replace `v'_m9=100*((5*`v' + 4*L.`v'   ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v' + 1*F4.`v')/19) if L2.`v'==.
 18.   replace `v'_m9=100*((5*`v'   ///
>                              + 4*F.`v' + 3*F2.`v' + 2*F3.`v' + 1*F4.`v')/15) if L.`v'==.
 19. end

. 
. cap program drop prepdata

. program define prepdata
  1.   syntax, indata(string) agerange(string) agegp(real)
  2.   use age educ4 unem empl yearmo wgt_composite if educ4==4 & `agerange' using `indata', clear
  3.   collapse (mean) empl unem [aw=wgt_composite], by(yearmo)
  4.   gen year=year(dofm(yearmo))
  5.   gen mon=month(dofm(yearmo))
  6.   seasadj empl
  7.   seasadj unem
  8.   smoother empl_sa
  9.   smoother unem_sa
 10.   gen noempl_sa_m7=100-empl_sa_m7
 11.   gen noempl_sa_m9=100-empl_sa_m9
 12.   gen agegp=`agegp'
 13.   gen recession=(yearmo>=ym(1981,7) & yearmo<=ym(1982,11)) | (yearmo>=ym(1990,7) & yearmo<=ym(1991,3)) | ///
>                 (yearmo>=ym(2001,3) & yearmo<=ym(2001,11)) | (yearmo>=ym(2007,12) & yearmo<=ym(2009,6))
 14. end

. 
. *** Prepare data
.  tempfile agegp0 agegp1 agegp2 agegp3 agegp4 agegp5

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=22 & age<=40") agegp(0)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        = 116.52
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.9202
                                                    Adj R-squared     = 0.9110
                                                    Root MSE          = 0.0046

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |    .000598   .0018891     0.32   0.752    -.0031148    .0043108
       1981  |   .0027381   .0018891     1.45   0.148    -.0009747     .006451
       1982  |  -.0018215   .0018891    -0.96   0.335    -.0055343    .0018914
       1983  |  -.0001538   .0018891    -0.08   0.935    -.0038666     .003559
       1984  |   .0119447   .0018891     6.32   0.000     .0082319    .0156575
       1985  |   .0198026   .0018891    10.48   0.000     .0160898    .0235154
       1986  |   .0228228   .0018891    12.08   0.000       .01911    .0265356
       1987  |   .0274019   .0018891    14.51   0.000      .023689    .0311147
       1988  |    .033139   .0018891    17.54   0.000     .0294262    .0368518
       1989  |   .0317073   .0018891    16.78   0.000     .0279945    .0354201
       1990  |    .028857   .0018891    15.28   0.000     .0251441    .0325698
       1991  |   .0178969   .0018891     9.47   0.000     .0141841    .0216097
       1992  |   .0183175   .0018891     9.70   0.000     .0146047    .0220303
       1993  |   .0227374   .0018891    12.04   0.000     .0190246    .0264502
       1994  |   .0262533   .0018891    13.90   0.000     .0225405    .0299661
       1995  |   .0248843   .0018891    13.17   0.000     .0211714    .0285971
       1996  |   .0250368   .0018891    13.25   0.000      .021324    .0287496
       1997  |   .0281292   .0018891    14.89   0.000     .0244164     .031842
       1998  |   .0244543   .0018891    12.94   0.000     .0207414    .0281671
       1999  |   .0229931   .0018891    12.17   0.000     .0192803    .0267059
       2000  |   .0203882   .0018891    10.79   0.000     .0166754     .024101
       2001  |    .007906   .0018891     4.19   0.000     .0041932    .0116188
       2002  |   .0008616   .0018891     0.46   0.649    -.0028512    .0045744
       2003  |  -.0085099   .0018891    -4.50   0.000    -.0122227   -.0047971
       2004  |  -.0008193   .0018891    -0.43   0.665    -.0045321    .0028935
       2005  |   .0026303   .0018891     1.39   0.165    -.0010825    .0063431
       2006  |   .0086594   .0018891     4.58   0.000     .0049466    .0123722
       2007  |   .0084592   .0018891     4.48   0.000     .0047464     .012172
       2008  |   .0064428   .0018891     3.41   0.001       .00273    .0101556
       2009  |  -.0120015   .0018891    -6.35   0.000    -.0157143   -.0082887
       2010  |  -.0165306   .0018891    -8.75   0.000    -.0202435   -.0128178
       2011  |  -.0166682   .0018891    -8.82   0.000    -.0203811   -.0129554
       2012  |  -.0115569   .0018891    -6.12   0.000    -.0152697   -.0078441
       2013  |  -.0110149   .0018891    -5.83   0.000    -.0147277   -.0073021
       2014  |  -.0071885   .0018891    -3.81   0.000    -.0109013   -.0034757
       2015  |  -.0055338   .0018891    -2.93   0.004    -.0092466    -.001821
       2016  |  -.0007724   .0018891    -0.41   0.683    -.0044852    .0029405
       2017  |   .0052963   .0018891     2.80   0.005     .0015834    .0090091
       2018  |   .0090034   .0018891     4.77   0.000     .0052906    .0127162
       2019  |   .0125058   .0018891     6.62   0.000      .008793    .0162187
             |
       _cons |   .8462156   .0013358   633.49   0.000     .8435902    .8488409
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 37.575            Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        = 109.04
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.9194
                                                    Adj R-squared     = 0.9101
                                                    Root MSE          = 0.0027

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |   .0017207   .0010903     1.58   0.115    -.0004221    .0038635
       1981  |    .004718   .0010903     4.33   0.000     .0025752    .0068607
       1982  |   .0137323   .0010903    12.60   0.000     .0115895     .015875
       1983  |   .0136551   .0010903    12.52   0.000     .0115123    .0157978
       1984  |   .0047823   .0010903     4.39   0.000     .0026395    .0069251
       1985  |   .0010807   .0010903     0.99   0.322     -.001062    .0032235
       1986  |   .0004643   .0010903     0.43   0.670    -.0016785    .0026071
       1987  |  -.0007429   .0010903    -0.68   0.496    -.0028856    .0013999
       1988  |  -.0051964   .0010903    -4.77   0.000    -.0073392   -.0030537
       1989  |  -.0037077   .0010903    -3.40   0.001    -.0058505   -.0015649
       1990  |  -.0020593   .0010903    -1.89   0.060    -.0042021    .0000835
       1991  |   .0065652   .0010903     6.02   0.000     .0044224     .008708
       1992  |   .0090536   .0010903     8.30   0.000     .0069108    .0111963
       1993  |   .0041852   .0010903     3.84   0.000     .0020425     .006328
       1994  |   .0004806   .0010903     0.44   0.660    -.0016622    .0026233
       1995  |  -.0005733   .0010903    -0.53   0.599    -.0027161    .0015695
       1996  |  -.0025361   .0010903    -2.33   0.020    -.0046789   -.0003933
       1997  |  -.0071319   .0010903    -6.54   0.000    -.0092746   -.0049891
       1998  |  -.0091466   .0010903    -8.39   0.000    -.0112894   -.0070038
       1999  |  -.0077808   .0010903    -7.14   0.000    -.0099236    -.005638
       2000  |  -.0095487   .0010903    -8.76   0.000    -.0116915    -.007406
       2001  |  -.0013606   .0010903    -1.25   0.213    -.0035034    .0007822
       2002  |   .0037458   .0010903     3.44   0.001      .001603    .0058886
       2003  |   .0063545   .0010903     5.83   0.000     .0042117    .0084972
       2004  |   .0017763   .0010903     1.63   0.104    -.0003665    .0039191
       2005  |  -.0023334   .0010903    -2.14   0.033    -.0044762   -.0001907
       2006  |  -.0053204   .0010903    -4.88   0.000    -.0074631   -.0031776
       2007  |  -.0042779   .0010903    -3.92   0.000    -.0064207   -.0021352
       2008  |   .0005713   .0010903     0.52   0.601    -.0015715    .0027141
       2009  |   .0226876   .0010903    20.81   0.000     .0205449    .0248304
       2010  |    .022439   .0010903    20.58   0.000     .0202963    .0245818
       2011  |   .0178522   .0010903    16.37   0.000     .0157094    .0199949
       2012  |   .0140587   .0010903    12.89   0.000     .0119159    .0162015
       2013  |   .0119071   .0010903    10.92   0.000     .0097643    .0140499
       2014  |     .00655   .0010903     6.01   0.000     .0044073    .0086928
       2015  |    .000077   .0010903     0.07   0.944    -.0020658    .0022197
       2016  |  -.0016135   .0010903    -1.48   0.140    -.0037562    .0005293
       2017  |  -.0031043   .0010903    -2.85   0.005    -.0052471   -.0009615
       2018  |  -.0047713   .0010903    -4.38   0.000    -.0069141   -.0026285
       2019  |   -.005203   .0010903    -4.77   0.000    -.0073457   -.0030602
             |
       _cons |   .0282321   .0007709    36.62   0.000     .0267169    .0297472
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 59.974            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp0'
file /tmp/St2868349.000004 saved as .dta format

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=22 & age<=30") agegp(1)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  79.52
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8841
                                                    Adj R-squared     = 0.8706
                                                    Root MSE          = 0.0070

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |  -.0035269   .0028704    -1.23   0.220    -.0091683    .0021145
       1981  |  -.0034799   .0028704    -1.21   0.226    -.0091213    .0021615
       1982  |  -.0116814   .0028704    -4.07   0.000    -.0173228     -.00604
       1983  |  -.0080066   .0028704    -2.79   0.006     -.013648   -.0023652
       1984  |   .0040343   .0028704     1.41   0.161    -.0016071    .0096757
       1985  |     .01747   .0028704     6.09   0.000     .0118286    .0231114
       1986  |   .0210651   .0028704     7.34   0.000     .0154237    .0267065
       1987  |   .0227947   .0028704     7.94   0.000     .0171533    .0284361
       1988  |   .0309083   .0028704    10.77   0.000     .0252669    .0365497
       1989  |    .028093   .0028704     9.79   0.000     .0224516    .0337344
       1990  |   .0227277   .0028704     7.92   0.000     .0170863    .0283691
       1991  |    .011147   .0028704     3.88   0.000     .0055056    .0167884
       1992  |   .0144709   .0028704     5.04   0.000     .0088295    .0201123
       1993  |   .0194557   .0028704     6.78   0.000     .0138143    .0250971
       1994  |   .0206546   .0028704     7.20   0.000     .0150131     .026296
       1995  |   .0231362   .0028704     8.06   0.000     .0174948    .0287777
       1996  |   .0237149   .0028704     8.26   0.000     .0180735    .0293563
       1997  |   .0303441   .0028704    10.57   0.000     .0247027    .0359855
       1998  |   .0263092   .0028704     9.17   0.000     .0206678    .0319506
       1999  |   .0193902   .0028704     6.76   0.000     .0137488    .0250316
       2000  |   .0229468   .0028704     7.99   0.000     .0173053    .0285882
       2001  |   .0066392   .0028704     2.31   0.021     .0009978    .0122806
       2002  |   .0009063   .0028704     0.32   0.752    -.0047351    .0065477
       2003  |  -.0087555   .0028704    -3.05   0.002    -.0143969    -.003114
       2004  |  -.0067042   .0028704    -2.34   0.020    -.0123456   -.0010628
       2005  |  -.0027903   .0028704    -0.97   0.332    -.0084318    .0028511
       2006  |   .0075498   .0028704     2.63   0.009     .0019084    .0131912
       2007  |   .0104711   .0028704     3.65   0.000     .0048297    .0161125
       2008  |    .001138   .0028704     0.40   0.692    -.0045034    .0067795
       2009  |  -.0231422   .0028704    -8.06   0.000    -.0287836   -.0175008
       2010  |  -.0283959   .0028704    -9.89   0.000    -.0340373   -.0227545
       2011  |  -.0283056   .0028704    -9.86   0.000     -.033947   -.0226642
       2012  |  -.0247923   .0028704    -8.64   0.000    -.0304337   -.0191509
       2013  |  -.0295441   .0028704   -10.29   0.000    -.0351855   -.0239027
       2014  |  -.0217431   .0028704    -7.57   0.000    -.0273845   -.0161017
       2015  |  -.0204935   .0028704    -7.14   0.000    -.0261349    -.014852
       2016  |  -.0140249   .0028704    -4.89   0.000    -.0196663   -.0083835
       2017  |  -.0052252   .0028704    -1.82   0.069    -.0108666    .0004162
       2018  |  -.0001151   .0028704    -0.04   0.968    -.0057565    .0055263
       2019  |  -.0000763   .0028704    -0.03   0.979    -.0057177    .0055651
             |
       _cons |    .842287   .0020297   414.98   0.000     .8382979    .8462761
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 15.920            Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  65.26
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8824
                                                    Adj R-squared     = 0.8688
                                                    Root MSE          = 0.0041

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |    .004632   .0016549     2.80   0.005     .0013794    .0078845
       1981  |   .0080408   .0016549     4.86   0.000     .0047882    .0112933
       1982  |   .0187151   .0016549    11.31   0.000     .0154625    .0219676
       1983  |   .0156615   .0016549     9.46   0.000      .012409    .0189141
       1984  |    .006699   .0016549     4.05   0.000     .0034465    .0099516
       1985  |    .001044   .0016549     0.63   0.528    -.0022085    .0042966
       1986  |   .0010021   .0016549     0.61   0.545    -.0022504    .0042547
       1987  |   .0007017   .0016549     0.42   0.672    -.0025509    .0039542
       1988  |  -.0061063   .0016549    -3.69   0.000    -.0093588   -.0028537
       1989  |   -.003888   .0016549    -2.35   0.019    -.0071405   -.0006354
       1990  |  -.0018769   .0016549    -1.13   0.257    -.0051294    .0013757
       1991  |   .0088098   .0016549     5.32   0.000     .0055572    .0120623
       1992  |    .008486   .0016549     5.13   0.000     .0052335    .0117386
       1993  |   .0037666   .0016549     2.28   0.023      .000514    .0070192
       1994  |   .0014598   .0016549     0.88   0.378    -.0017928    .0047124
       1995  |  -.0007111   .0016549    -0.43   0.668    -.0039637    .0025414
       1996  |   -.001957   .0016549    -1.18   0.238    -.0052095    .0012956
       1997  |  -.0097385   .0016549    -5.88   0.000    -.0129911    -.006486
       1998  |  -.0109837   .0016549    -6.64   0.000    -.0142363   -.0077311
       1999  |  -.0072029   .0016549    -4.35   0.000    -.0104555   -.0039504
       2000  |  -.0119129   .0016549    -7.20   0.000    -.0151654   -.0086603
       2001  |  -.0008515   .0016549    -0.51   0.607    -.0041041     .002401
       2002  |   .0033831   .0016549     2.04   0.042     .0001306    .0066357
       2003  |   .0046101   .0016549     2.79   0.006     .0013575    .0078626
       2004  |   .0025612   .0016549     1.55   0.122    -.0006914    .0058137
       2005  |  -.0013327   .0016549    -0.81   0.421    -.0045853    .0019198
       2006  |    -.00565   .0016549    -3.41   0.001    -.0089025   -.0023974
       2007  |  -.0043193   .0016549    -2.61   0.009    -.0075719   -.0010667
       2008  |   .0011401   .0016549     0.69   0.491    -.0021125    .0043927
       2009  |   .0239941   .0016549    14.50   0.000     .0207415    .0272467
       2010  |   .0284179   .0016549    17.17   0.000     .0251653    .0316704
       2011  |    .021298   .0016549    12.87   0.000     .0180454    .0245506
       2012  |   .0162203   .0016549     9.80   0.000     .0129678    .0194729
       2013  |   .0169422   .0016549    10.24   0.000     .0136896    .0201947
       2014  |   .0092272   .0016549     5.58   0.000     .0059746    .0124798
       2015  |   .0016465   .0016549     0.99   0.320    -.0016061    .0048991
       2016  |  -.0012077   .0016549    -0.73   0.466    -.0044603    .0020449
       2017  |  -.0024859   .0016549    -1.50   0.134    -.0057384    .0007667
       2018  |  -.0047044   .0016549    -2.84   0.005     -.007957   -.0014519
       2019  |  -.0043498   .0016549    -2.63   0.009    -.0076024   -.0010973
             |
       _cons |   .0344258   .0011702    29.42   0.000     .0321259    .0367257
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 62.773            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp1'
file /tmp/St2868349.000005 saved as .dta format

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=31 & age<=40") agegp(2)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  61.38
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8656
                                                    Adj R-squared     = 0.8500
                                                    Root MSE          = 0.0054

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |   .0044938   .0021852     2.06   0.040     .0001991    .0087884
       1981  |   .0079366   .0021852     3.63   0.000      .003642    .0122313
       1982  |   .0061501   .0021852     2.81   0.005     .0018554    .0104447
       1983  |   .0056948   .0021852     2.61   0.009     .0014002    .0099895
       1984  |   .0173001   .0021852     7.92   0.000     .0130054    .0215947
       1985  |   .0204922   .0021852     9.38   0.000     .0161975    .0247868
       1986  |   .0228258   .0021852    10.45   0.000     .0185312    .0271205
       1987  |   .0293848   .0021852    13.45   0.000     .0250902    .0336795
       1988  |   .0333253   .0021852    15.25   0.000     .0290307      .03762
       1989  |   .0328116   .0021852    15.02   0.000      .028517    .0371062
       1990  |   .0318335   .0021852    14.57   0.000     .0275389    .0361281
       1991  |   .0213855   .0021852     9.79   0.000     .0170908    .0256801
       1992  |   .0195053   .0021852     8.93   0.000     .0152107       .0238
       1993  |   .0235095   .0021852    10.76   0.000     .0192149    .0278042
       1994  |   .0286866   .0021852    13.13   0.000     .0243919    .0329812
       1995  |   .0245509   .0021852    11.24   0.000     .0202563    .0288456
       1996  |   .0244358   .0021852    11.18   0.000     .0201412    .0287305
       1997  |   .0250172   .0021852    11.45   0.000     .0207226    .0293118
       1998  |   .0215821   .0021852     9.88   0.000     .0172875    .0258768
       1999  |   .0239849   .0021852    10.98   0.000     .0196902    .0282795
       2000  |     .01704   .0021852     7.80   0.000     .0127453    .0213346
       2001  |   .0071123   .0021852     3.25   0.001     .0028177    .0114069
       2002  |   -.000869   .0021852    -0.40   0.691    -.0051637    .0034256
       2003  |  -.0099678   .0021852    -4.56   0.000    -.0142625   -.0056732
       2004  |   .0017604   .0021852     0.81   0.421    -.0025342     .006055
       2005  |   .0049936   .0021852     2.29   0.023      .000699    .0092883
       2006  |   .0078973   .0021852     3.61   0.000     .0036027    .0121919
       2007  |   .0055873   .0021852     2.56   0.011     .0012927     .009882
       2008  |   .0091984   .0021852     4.21   0.000     .0049038    .0134931
       2009  |  -.0045751   .0021852    -2.09   0.037    -.0088697   -.0002804
       2010  |  -.0081934   .0021852    -3.75   0.000     -.012488   -.0038987
       2011  |  -.0083426   .0021852    -3.82   0.000    -.0126372   -.0040479
       2012  |  -.0018421   .0021852    -0.84   0.400    -.0061367    .0024526
       2013  |   .0027724   .0021852     1.27   0.205    -.0015223     .007067
       2014  |    .003443   .0021852     1.58   0.116    -.0008516    .0077377
       2015  |    .005609   .0021852     2.57   0.011     .0013143    .0099036
       2016  |   .0086745   .0021852     3.97   0.000     .0043798    .0129691
       2017  |    .012574   .0021852     5.75   0.000     .0082793    .0168686
       2018  |   .0150444   .0021852     6.88   0.000     .0107498    .0193391
       2019  |    .021174   .0021852     9.69   0.000     .0168793    .0254686
             |
       _cons |    .850553   .0015451   550.47   0.000     .8475162    .8535898
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 34.361            Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  81.01
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8873
                                                    Adj R-squared     = 0.8743
                                                    Root MSE          = 0.0027

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |  -.0008862   .0011028    -0.80   0.422    -.0030536    .0012813
       1981  |   .0025842   .0011028     2.34   0.020     .0004167    .0047517
       1982  |    .010595   .0011028     9.61   0.000     .0084276    .0127625
       1983  |    .013378   .0011028    12.13   0.000     .0112105    .0155455
       1984  |    .004918   .0011028     4.46   0.000     .0027505    .0070855
       1985  |   .0028769   .0011028     2.61   0.009     .0007094    .0050444
       1986  |   .0021069   .0011028     1.91   0.057    -.0000606    .0042744
       1987  |   .0003606   .0011028     0.33   0.744    -.0018069     .002528
       1988  |  -.0023362   .0011028    -2.12   0.035    -.0045037   -.0001687
       1989  |  -.0012799   .0011028    -1.16   0.246    -.0034474    .0008876
       1990  |    .000104   .0011028     0.09   0.925    -.0020635    .0022715
       1991  |   .0071957   .0011028     6.52   0.000     .0050282    .0093632
       1992  |   .0118607   .0011028    10.75   0.000     .0096933    .0140282
       1993  |   .0068916   .0011028     6.25   0.000     .0047241    .0090591
       1994  |    .002174   .0011028     1.97   0.049     6.49e-06    .0043415
       1995  |   .0019354   .0011028     1.75   0.080    -.0002321    .0041029
       1996  |  -.0005953   .0011028    -0.54   0.590    -.0027628    .0015722
       1997  |  -.0029245   .0011028    -2.65   0.008     -.005092    -.000757
       1998  |  -.0054352   .0011028    -4.93   0.000    -.0076027   -.0032677
       1999  |  -.0057676   .0011028    -5.23   0.000     -.007935   -.0036001
       2000  |  -.0055277   .0011028    -5.01   0.000    -.0076952   -.0033602
       2001  |   .0008148   .0011028     0.74   0.460    -.0013527    .0029823
       2002  |   .0065585   .0011028     5.95   0.000      .004391     .008726
       2003  |   .0100593   .0011028     9.12   0.000     .0078918    .0122268
       2004  |   .0036685   .0011028     3.33   0.001      .001501     .005836
       2005  |  -.0007159   .0011028    -0.65   0.517    -.0028834    .0014515
       2006  |  -.0028526   .0011028    -2.59   0.010    -.0050201   -.0006852
       2007  |  -.0022359   .0011028    -2.03   0.043    -.0044034   -.0000684
       2008  |   .0021985   .0011028     1.99   0.047      .000031     .004366
       2009  |   .0237257   .0011028    21.51   0.000     .0215582    .0258932
       2010  |   .0195492   .0011028    17.73   0.000     .0173817    .0217166
       2011  |   .0168932   .0011028    15.32   0.000     .0147257    .0190607
       2012  |    .014132   .0011028    12.81   0.000     .0119645    .0162995
       2013  |   .0098519   .0011028     8.93   0.000     .0076844    .0120194
       2014  |    .006293   .0011028     5.71   0.000     .0041255    .0084605
       2015  |   .0006508   .0011028     0.59   0.555    -.0015167    .0028183
       2016  |  -.0000311   .0011028    -0.03   0.978    -.0021986    .0021364
       2017  |  -.0017325   .0011028    -1.57   0.117       -.0039    .0004349
       2018  |  -.0029028   .0011028    -2.63   0.009    -.0050703   -.0007353
       2019  |  -.0039004   .0011028    -3.54   0.000    -.0060679   -.0017329
             |
       _cons |   .0214611   .0007798    27.52   0.000     .0199285    .0229937
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 20.426            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp2'
file /tmp/St2868349.000006 saved as .dta format

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=22 & age<=25") agegp(3)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  46.55
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8239
                                                    Adj R-squared     = 0.8035
                                                    Root MSE          = 0.0118

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |  -.0095855   .0048352    -1.98   0.048    -.0190885   -.0000826
       1981  |  -.0113624   .0048352    -2.35   0.019    -.0208653   -.0018595
       1982  |  -.0202201   .0048352    -4.18   0.000     -.029723   -.0107171
       1983  |  -.0259264   .0048352    -5.36   0.000    -.0354293   -.0164235
       1984  |  -.0116078   .0048352    -2.40   0.017    -.0211107   -.0021049
       1985  |   .0043893   .0048352     0.91   0.364    -.0051136    .0138923
       1986  |   .0109095   .0048352     2.26   0.025     .0014066    .0204124
       1987  |   .0173628   .0048352     3.59   0.000     .0078599    .0268658
       1988  |    .030772   .0048352     6.36   0.000      .021269    .0402749
       1989  |   .0165832   .0048352     3.43   0.001     .0070803    .0260861
       1990  |   .0104878   .0048352     2.17   0.031     .0009849    .0199907
       1991  |  -.0089518   .0048352    -1.85   0.065    -.0184547    .0005512
       1992  |    .006898   .0048352     1.43   0.154    -.0026049    .0164009
       1993  |   .0079671   .0048352     1.65   0.100    -.0015358      .01747
       1994  |   .0078431   .0048352     1.62   0.106    -.0016598     .017346
       1995  |   .0017399   .0048352     0.36   0.719     -.007763    .0112428
       1996  |   .0036654   .0048352     0.76   0.449    -.0058375    .0131683
       1997  |   .0192711   .0048352     3.99   0.000     .0097681     .028774
       1998  |   .0086048   .0048352     1.78   0.076    -.0008982    .0181077
       1999  |   .0023855   .0048352     0.49   0.622    -.0071174    .0118884
       2000  |   .0074608   .0048352     1.54   0.124    -.0020421    .0169637
       2001  |  -.0051601   .0048352    -1.07   0.286     -.014663    .0043428
       2002  |  -.0168519   .0048352    -3.49   0.001    -.0263548    -.007349
       2003  |  -.0266309   .0048352    -5.51   0.000    -.0361338   -.0171279
       2004  |   -.022814   .0048352    -4.72   0.000    -.0323169   -.0133111
       2005  |  -.0144414   .0048352    -2.99   0.003    -.0239443   -.0049385
       2006  |  -.0074992   .0048352    -1.55   0.122    -.0170021    .0020037
       2007  |  -.0031907   .0048352    -0.66   0.510    -.0126936    .0063123
       2008  |  -.0256152   .0048352    -5.30   0.000    -.0351181   -.0161123
       2009  |  -.0450034   .0048352    -9.31   0.000    -.0545063   -.0355005
       2010  |  -.0485649   .0048352   -10.04   0.000    -.0580678    -.039062
       2011  |  -.0532204   .0048352   -11.01   0.000    -.0627233   -.0437175
       2012  |   -.053428   .0048352   -11.05   0.000     -.062931   -.0439251
       2013  |   -.062272   .0048352   -12.88   0.000    -.0717749   -.0527691
       2014  |  -.0410926   .0048352    -8.50   0.000    -.0505955   -.0315897
       2015  |  -.0468417   .0048352    -9.69   0.000    -.0563447   -.0373388
       2016  |  -.0408134   .0048352    -8.44   0.000    -.0503164   -.0313105
       2017  |   -.032935   .0048352    -6.81   0.000     -.042438   -.0234321
       2018  |  -.0309392   .0048352    -6.40   0.000    -.0404421   -.0214363
       2019  |  -.0283252   .0048352    -5.86   0.000    -.0378281   -.0188223
             |
       _cons |    .825422    .003419   241.42   0.000     .8187024    .8321416
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 17.837            Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  31.14
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.8227
                                                    Adj R-squared     = 0.8022
                                                    Root MSE          = 0.0072

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |   .0072519   .0029305     2.47   0.014     .0014923    .0130114
       1981  |    .013471   .0029305     4.60   0.000     .0077114    .0192306
       1982  |   .0245577   .0029305     8.38   0.000     .0187981    .0303173
       1983  |   .0243413   .0029305     8.31   0.000     .0185817    .0301009
       1984  |    .008827   .0029305     3.01   0.003     .0030674    .0145866
       1985  |    .002851   .0029305     0.97   0.331    -.0029086    .0086106
       1986  |   .0053514   .0029305     1.83   0.069    -.0004082     .011111
       1987  |    .002491   .0029305     0.85   0.396    -.0032685    .0082506
       1988  |  -.0052338   .0029305    -1.79   0.075    -.0109933    .0005258
       1989  |  -.0032841   .0029305    -1.12   0.263    -.0090437    .0024755
       1990  |  -.0009368   .0029305    -0.32   0.749    -.0066963    .0048228
       1991  |   .0139509   .0029305     4.76   0.000     .0081914    .0197105
       1992  |   .0090242   .0029305     3.08   0.002     .0032647    .0147838
       1993  |   .0056452   .0029305     1.93   0.055    -.0001144    .0114048
       1994  |   .0011806   .0029305     0.40   0.687     -.004579    .0069402
       1995  |   .0023156   .0029305     0.79   0.430     -.003444    .0080752
       1996  |  -.0002875   .0029305    -0.10   0.922     -.006047    .0054721
       1997  |  -.0137392   .0029305    -4.69   0.000    -.0194988   -.0079796
       1998  |  -.0105877   .0029305    -3.61   0.000    -.0163473   -.0048281
       1999  |  -.0044952   .0029305    -1.53   0.126    -.0102548    .0012644
       2000  |    -.01114   .0029305    -3.80   0.000    -.0168996   -.0053804
       2001  |   .0020464   .0029305     0.70   0.485    -.0037132    .0078059
       2002  |   .0049017   .0029305     1.67   0.095    -.0008579    .0106613
       2003  |   .0057803   .0029305     1.97   0.049     .0000207    .0115399
       2004  |   .0052798   .0029305     1.80   0.072    -.0004798    .0110394
       2005  |   .0002838   .0029305     0.10   0.923    -.0054758    .0060433
       2006  |  -.0064032   .0029305    -2.19   0.029    -.0121628   -.0006436
       2007  |  -.0024609   .0029305    -0.84   0.402    -.0082205    .0032986
       2008  |   .0045878   .0029305     1.57   0.118    -.0011718    .0103474
       2009  |   .0292662   .0029305     9.99   0.000     .0235066    .0350258
       2010  |   .0340347   .0029305    11.61   0.000     .0282751    .0397943
       2011  |   .0317838   .0029305    10.85   0.000     .0260242    .0375434
       2012  |   .0212195   .0029305     7.24   0.000       .01546    .0269791
       2013  |   .0236286   .0029305     8.06   0.000      .017869    .0293882
       2014  |   .0169103   .0029305     5.77   0.000     .0111507    .0226699
       2015  |   .0064073   .0029305     2.19   0.029     .0006477    .0121668
       2016  |   .0030946   .0029305     1.06   0.292     -.002665    .0088542
       2017  |   .0009825   .0029305     0.34   0.738    -.0047771     .006742
       2018  |  -.0004587   .0029305    -0.16   0.876    -.0062183    .0053009
       2019  |  -.0017876   .0029305    -0.61   0.542    -.0075472     .003972
             |
       _cons |   .0450908   .0020722    21.76   0.000     .0410182    .0491635
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 72.369            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp3'
file /tmp/St2868349.000007 saved as .dta format

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=26 & age<=29") agegp(4)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  38.51
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.7840
                                                    Adj R-squared     = 0.7590
                                                    Root MSE          = 0.0091

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |  -.0010901    .003699    -0.29   0.768      -.00836    .0061798
       1981  |   .0020301    .003699     0.55   0.583    -.0052398       .0093
       1982  |  -.0073291    .003699    -1.98   0.048     -.014599   -.0000592
       1983  |   .0002294    .003699     0.06   0.951    -.0070405    .0074993
       1984  |   .0106201    .003699     2.87   0.004     .0033502      .01789
       1985  |   .0255757    .003699     6.91   0.000     .0183058    .0328456
       1986  |    .027438    .003699     7.42   0.000     .0201681    .0347079
       1987  |    .022871    .003699     6.18   0.000     .0156011    .0301409
       1988  |   .0290133    .003699     7.84   0.000     .0217434    .0362832
       1989  |   .0327105    .003699     8.84   0.000     .0254406    .0399805
       1990  |   .0305206    .003699     8.25   0.000     .0232507    .0377905
       1991  |   .0209468    .003699     5.66   0.000     .0136769    .0282167
       1992  |   .0170136    .003699     4.60   0.000     .0097437    .0242835
       1993  |   .0249467    .003699     6.74   0.000     .0176767    .0322166
       1994  |    .028743    .003699     7.77   0.000     .0214731    .0360129
       1995  |    .032246    .003699     8.72   0.000     .0249761    .0395159
       1996  |   .0336128    .003699     9.09   0.000     .0263429    .0408827
       1997  |   .0326456    .003699     8.83   0.000     .0253757    .0399155
       1998  |   .0298857    .003699     8.08   0.000     .0226158    .0371556
       1999  |   .0250343    .003699     6.77   0.000     .0177644    .0323042
       2000  |   .0316959    .003699     8.57   0.000      .024426    .0389658
       2001  |   .0113789    .003699     3.08   0.002      .004109    .0186488
       2002  |    .010501    .003699     2.84   0.005      .003231    .0177709
       2003  |   .0013685    .003699     0.37   0.712    -.0059015    .0086384
       2004  |   .0025223    .003699     0.68   0.496    -.0047476    .0097922
       2005  |   .0040556    .003699     1.10   0.274    -.0032143    .0113255
       2006  |   .0134366    .003699     3.63   0.000     .0061667    .0207065
       2007  |   .0166762    .003699     4.51   0.000     .0094063    .0239461
       2008  |   .0167294    .003699     4.52   0.000     .0094595    .0239993
       2009  |  -.0091101    .003699    -2.46   0.014      -.01638   -.0018402
       2010  |    -.01743    .003699    -4.71   0.000    -.0246999   -.0101601
       2011  |  -.0174353    .003699    -4.71   0.000    -.0247053   -.0101654
       2012  |  -.0128875    .003699    -3.48   0.001    -.0201574   -.0056176
       2013  |  -.0117071    .003699    -3.16   0.002    -.0189771   -.0044372
       2014  |  -.0178898    .003699    -4.84   0.000    -.0251597   -.0106199
       2015  |  -.0103985    .003699    -2.81   0.005    -.0176684   -.0031286
       2016  |  -.0026909    .003699    -0.73   0.467    -.0099608     .004579
       2017  |   .0088716    .003699     2.40   0.017     .0016017    .0161415
       2018  |   .0130837    .003699     3.54   0.000     .0058138    .0203536
       2019  |   .0126004    .003699     3.41   0.001     .0053305    .0198703
             |
       _cons |   .8539741   .0026156   326.49   0.000     .8488335    .8591147
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 5.176             Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        =  38.61
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.7925
                                                    Adj R-squared     = 0.7685
                                                    Root MSE          = 0.0046

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |   .0031124   .0018619     1.67   0.095    -.0005468    .0067717
       1981  |   .0044736   .0018619     2.40   0.017     .0008143    .0081328
       1982  |   .0153622   .0018619     8.25   0.000     .0117029    .0190214
       1983  |   .0106924   .0018619     5.74   0.000     .0070332    .0143517
       1984  |   .0050749   .0018619     2.73   0.007     .0014156    .0087341
       1985  |  -.0006778   .0018619    -0.36   0.716    -.0043371    .0029814
       1986  |  -.0008697   .0018619    -0.47   0.641     -.004529    .0027895
       1987  |   .0003118   .0018619     0.17   0.867    -.0033475     .003971
       1988  |  -.0063832   .0018619    -3.43   0.001    -.0100425    -.002724
       1989  |  -.0039096   .0018619    -2.10   0.036    -.0075689   -.0002504
       1990  |  -.0032266   .0018619    -1.73   0.084    -.0068858    .0004327
       1991  |   .0063996   .0018619     3.44   0.001     .0027404    .0100589
       1992  |   .0086613   .0018619     4.65   0.000      .005002    .0123205
       1993  |   .0032818   .0018619     1.76   0.079    -.0003775     .006941
       1994  |   .0023845   .0018619     1.28   0.201    -.0012748    .0060437
       1995  |  -.0015962   .0018619    -0.86   0.392    -.0052555     .002063
       1996  |  -.0023536   .0018619    -1.26   0.207    -.0060128    .0013057
       1997  |  -.0070491   .0018619    -3.79   0.000    -.0107084   -.0033899
       1998  |   -.009696   .0018619    -5.21   0.000    -.0133553   -.0060368
       1999  |  -.0077024   .0018619    -4.14   0.000    -.0113617   -.0040432
       2000  |  -.0124684   .0018619    -6.70   0.000    -.0161276   -.0088091
       2001  |  -.0026556   .0018619    -1.43   0.154    -.0063148    .0010037
       2002  |   .0035438   .0018619     1.90   0.058    -.0001155     .007203
       2003  |   .0042747   .0018619     2.30   0.022     .0006154    .0079339
       2004  |   .0004225   .0018619     0.23   0.821    -.0032368    .0040817
       2005  |  -.0014967   .0018619    -0.80   0.422    -.0051559    .0021626
       2006  |  -.0050754   .0018619    -2.73   0.007    -.0087347   -.0014162
       2007  |   -.004515   .0018619    -2.42   0.016    -.0081742   -.0008557
       2008  |  -.0017802   .0018619    -0.96   0.340    -.0054394    .0018791
       2009  |   .0209987   .0018619    11.28   0.000     .0173394    .0246579
       2010  |    .023932   .0018619    12.85   0.000     .0202727    .0275912
       2011  |   .0170553   .0018619     9.16   0.000     .0133961    .0207146
       2012  |   .0130447   .0018619     7.01   0.000     .0093855     .016704
       2013  |   .0129395   .0018619     6.95   0.000     .0092803    .0165988
       2014  |   .0049289   .0018619     2.65   0.008     .0012696    .0085881
       2015  |  -.0000458   .0018619    -0.02   0.980    -.0037051    .0036134
       2016  |  -.0025836   .0018619    -1.39   0.166    -.0062428    .0010757
       2017  |  -.0044626   .0018619    -2.40   0.017    -.0081219   -.0008034
       2018  |   -.005512   .0018619    -2.96   0.003    -.0091713   -.0018528
       2019  |  -.0044081   .0018619    -2.37   0.018    -.0080673   -.0007488
             |
       _cons |   .0289919   .0013165    22.02   0.000     .0264045    .0315794
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 12.386            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp4'
file /tmp/St2868349.000008 saved as .dta format

.  prepdata, indata("`scratch'/extractcps.dta") agerange("age>=40 & age<=55") agegp(4)

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        = 174.09
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.9455
                                                    Adj R-squared     = 0.9392
                                                    Root MSE          = 0.0041

------------------------------------------------------------------------------
        empl | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |   .0079176   .0016623     4.76   0.000     .0046505    .0111846
       1981  |   .0118374   .0016623     7.12   0.000     .0085704    .0151045
       1982  |   .0079677   .0016623     4.79   0.000     .0047006    .0112348
       1983  |   .0047877   .0016623     2.88   0.004     .0015206    .0080547
       1984  |   .0156121   .0016623     9.39   0.000      .012345    .0188792
       1985  |   .0230581   .0016623    13.87   0.000     .0197911    .0263252
       1986  |    .025123   .0016623    15.11   0.000     .0218559      .02839
       1987  |   .0276113   .0016623    16.61   0.000     .0243442    .0308783
       1988  |   .0324135   .0016623    19.50   0.000     .0291464    .0356805
       1989  |   .0320424   .0016623    19.28   0.000     .0287754    .0353095
       1990  |   .0338537   .0016623    20.37   0.000     .0305866    .0371208
       1991  |   .0294149   .0016623    17.70   0.000     .0261478    .0326819
       1992  |   .0261563   .0016623    15.73   0.000     .0228892    .0294234
       1993  |   .0278118   .0016623    16.73   0.000     .0245447    .0310789
       1994  |   .0295951   .0016623    17.80   0.000      .026328    .0328621
       1995  |   .0272895   .0016623    16.42   0.000     .0240224    .0305565
       1996  |   .0276882   .0016623    16.66   0.000     .0244211    .0309552
       1997  |   .0301387   .0016623    18.13   0.000     .0268716    .0334058
       1998  |   .0275359   .0016623    16.56   0.000     .0242688    .0308029
       1999  |   .0253558   .0016623    15.25   0.000     .0220887    .0286228
       2000  |    .023275   .0016623    14.00   0.000     .0200079    .0265421
       2001  |     .01534   .0016623     9.23   0.000     .0120729     .018607
       2002  |   .0049222   .0016623     2.96   0.003     .0016552    .0081893
       2003  |   .0025679   .0016623     1.54   0.123    -.0006991     .005835
       2004  |   .0003637   .0016623     0.22   0.827    -.0029033    .0036308
       2005  |    .006544   .0016623     3.94   0.000     .0032769    .0098111
       2006  |   .0097058   .0016623     5.84   0.000     .0064387    .0129728
       2007  |   .0104587   .0016623     6.29   0.000     .0071916    .0137258
       2008  |   .0065328   .0016623     3.93   0.000     .0032658    .0097999
       2009  |  -.0126904   .0016623    -7.63   0.000    -.0159574   -.0094233
       2010  |  -.0190674   .0016623   -11.47   0.000    -.0223345   -.0158003
       2011  |  -.0157886   .0016623    -9.50   0.000    -.0190557   -.0125215
       2012  |  -.0125762   .0016623    -7.57   0.000    -.0158432   -.0093091
       2013  |  -.0143814   .0016623    -8.65   0.000    -.0176484   -.0111143
       2014  |  -.0075745   .0016623    -4.56   0.000    -.0108416   -.0043074
       2015  |  -.0033789   .0016623    -2.03   0.043     -.006646   -.0001119
       2016  |  -.0015562   .0016623    -0.94   0.350    -.0048233    .0017109
       2017  |  -.0019091   .0016623    -1.15   0.251    -.0051761     .001358
       2018  |   .0025228   .0016623     1.52   0.130    -.0007443    .0057898
       2019  |   .0049376   .0016623     2.97   0.003     .0016705    .0082047
             |
       _cons |   .8620927   .0011754   733.42   0.000     .8597826    .8644029
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 60.779            Prob > F = 0.000

Linear regression, absorbing indicators             Number of obs     =    492
Absorbed variable: mon                              No. of categories =     12
                                                    F(40, 440)        = 101.06
                                                    Prob > F          = 0.0000
                                                    R-squared         = 0.9053
                                                    Adj R-squared     = 0.8943
                                                    Root MSE          = 0.0025

------------------------------------------------------------------------------
        unem | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
        year |
       1980  |  -.0002898   .0010314    -0.28   0.779    -.0023169    .0017374
       1981  |   .0022726   .0010314     2.20   0.028     .0002454    .0042998
       1982  |   .0058963   .0010314     5.72   0.000     .0038691    .0079235
       1983  |   .0083272   .0010314     8.07   0.000     .0063001    .0103544
       1984  |   .0044271   .0010314     4.29   0.000     .0023999    .0064543
       1985  |   .0031332   .0010314     3.04   0.003      .001106    .0051604
       1986  |   .0038273   .0010314     3.71   0.000     .0018001    .0058545
       1987  |   .0036294   .0010314     3.52   0.000     .0016023    .0056566
       1988  |   .0011853   .0010314     1.15   0.251    -.0008419    .0032124
       1989  |   .0031396   .0010314     3.04   0.002     .0011124    .0051668
       1990  |   .0029267   .0010314     2.84   0.005     .0008995    .0049539
       1991  |   .0069129   .0010314     6.70   0.000     .0048857      .00894
       1992  |   .0122783   .0010314    11.90   0.000     .0102511    .0143055
       1993  |   .0106966   .0010314    10.37   0.000     .0086694    .0127237
       1994  |   .0076498   .0010314     7.42   0.000     .0056227     .009677
       1995  |    .006372   .0010314     6.18   0.000     .0043449    .0083992
       1996  |   .0045433   .0010314     4.40   0.000     .0025162    .0065705
       1997  |    .002857   .0010314     2.77   0.006     .0008299    .0048842
       1998  |    .000379   .0010314     0.37   0.713    -.0016482    .0024062
       1999  |  -.0004778   .0010314    -0.46   0.643     -.002505    .0015494
       2000  |  -.0007588   .0010314    -0.74   0.462     -.002786    .0012684
       2001  |   .0038883   .0010314     3.77   0.000     .0018612    .0059155
       2002  |   .0111405   .0010314    10.80   0.000     .0091134    .0131677
       2003  |   .0116297   .0010314    11.28   0.000     .0096025    .0136569
       2004  |   .0096168   .0010314     9.32   0.000     .0075896    .0116439
       2005  |   .0052009   .0010314     5.04   0.000     .0031737    .0072281
       2006  |   .0022815   .0010314     2.21   0.027     .0002543    .0043087
       2007  |   .0014334   .0010314     1.39   0.165    -.0005937    .0034606
       2008  |   .0066121   .0010314     6.41   0.000     .0045849    .0086392
       2009  |   .0272747   .0010314    26.44   0.000     .0252475    .0293018
       2010  |   .0292417   .0010314    28.35   0.000     .0272146    .0312689
       2011  |   .0242457   .0010314    23.51   0.000     .0222185    .0262729
       2012  |   .0205491   .0010314    19.92   0.000      .018522    .0225763
       2013  |   .0180302   .0010314    17.48   0.000      .016003    .0200574
       2014  |      .0116   .0010314    11.25   0.000     .0095728    .0136272
       2015  |   .0062333   .0010314     6.04   0.000     .0042062    .0082605
       2016  |   .0059311   .0010314     5.75   0.000     .0039039    .0079583
       2017  |   .0040323   .0010314     3.91   0.000     .0020052    .0060595
       2018  |   .0028851   .0010314     2.80   0.005     .0008579    .0049123
       2019  |   .0023904   .0010314     2.32   0.021     .0003632    .0044176
             |
       _cons |   .0170558   .0007293    23.39   0.000     .0156224    .0184892
------------------------------------------------------------------------------
F test of absorbed indicators: F(11, 440) = 14.930            Prob > F = 0.000

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

Time variable: yearmo, 1979m1 to 2019m12
        Delta: 1 month
(6 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(8 missing values generated)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)
(1 real change made)

.  save `agegp5'
file /tmp/St2868349.000009 saved as .dta format

.  use `agegp0'

.  append using `agegp1'

.  append using `agegp2'

.  append using `agegp3'

.  append using `agegp4'

.  append using `agegp5'

.  label def agegp_l 0 "22-40" 1 "22-30" 2 "31-40" 3 "22-25" 4 "26-29" 5 "40-55"

.  label values agegp agegp_l

.  
. * Figure 2 - 22-40 Graduate Non-Employment & UR National rate
. gen altrecession=8*recession+10 

. twoway area altrecession yearmo if agegp==0, yaxis(2) color(gs13) base(10) || ///
>  line  noempl_sa_m7 yearmo if agegp==0, yaxis(2) ytitle("Non-employment rate", axis(2)) ylabel(10(2)18, axis(2)) || ///
>  line  unem_sa_m7  yearmo if agegp==0, yaxis(3) ytitle("Unemployment rate", axis(3)) lpattern(dash)  || ///
>  , xlabel(246 "1980" 306 "1985" 366 "1990" 426 "1995" 486 "2000" 546 "2005" 606 "2010" 666 "2015" 726 "2020") ///
>  xtitle("")  title("Non-employment and unemployment, college graduates aged 22-40", size(medsmall)) legend(order(2  "Non-employmen
> t rate (l. axis)" 3 "Unemployment rate (r. axis)") size(small)) ///
> saving("`output'/`dofile'_base.gph", replace)
file /accounts/projects/jr_ra/GRscarring/erratum/results/fig_ur_age_base.gph saved

. 
. 
. 
. * Alternative -- non-employment for 25-30 and 31-40
. replace altrecession=10+10*recession
(324 real changes made)

. twoway area altrecession yearmo if agegp==1, yaxis(1) color(gs13) base(10) || ///
>  line  noempl_sa_m7 yearmo if agegp==1, yaxis(1) lstyle(p1) || ///
>  line noempl_sa_m7 yearmo if agegp==2, yaxis(1) lstyle(p2) lpattern(dash) || ///
>  , ytitle("Non-employment rate") yla(10 (2) 20) ///
>   xlabel(246 "1980" 306 "1985" 366 "1990" 426 "1995" 486 "2000" 546 "2005" 606 "2010" 666 "2015" 726 "2020") ///
>  xtitle("")  title("Non-employment rate, by age", size(medsmall)) legend(order(2 "Age 22-30" 3 "Age 31-40") col(2)) ///
> saving("`output'/`dofile'_agegp2_noempl.gph", replace)
file /accounts/projects/jr_ra/GRscarring/erratum/results/fig_ur_age_agegp2_noempl.gph saved

. 
. gen altrecessionB=7*recession

. twoway area altrecessionB yearmo if agegp==1, yaxis(1) color(gs13) || ///
>  line  unem_sa_m7 yearmo if agegp==1, yaxis(1) lstyle(p1) || ///
>  line unem_sa_m7 yearmo if agegp==2, yaxis(1) lstyle(p2)  lpattern(dash) || ///
>  , ytitle("Non-employment rate") yla(0 (1) 6) ///
>   xlabel(246 "1980" 306 "1985" 366 "1990" 426 "1995" 486 "2000" 546 "2005" 606 "2010" 666 "2015" 726 "2020") ///
>  xtitle("")  title("Unemployment rate, by age", size(medsmall)) legend(order(2 "Age 22-30" 3 "Age 31-40") col(2)) ///
> saving("`output'/`dofile'_agegp2_unem.gph", replace)
file /accounts/projects/jr_ra/GRscarring/erratum/results/fig_ur_age_agegp2_unem.gph saved

. 
. gen altrecessionC=12*recession+78

. twoway area altrecessionC yearmo if agegp==1, base(78) yaxis(1) color(gs13) || ///
>  line  empl_sa_m7 yearmo if agegp==1, yaxis(1) lstyle(p1) || ///
>  line empl_sa_m7 yearmo if agegp==2, yaxis(1) lstyle(p2)  lpattern(dash) || ///
>  , ytitle("Employment rate") yla(78 (2) 90) ///
>   xlabel(246 "1980" 306 "1985" 366 "1990" 426 "1995" 486 "2000" 546 "2005" 606 "2010" 666 "2015" 726 "2020") ///
>  xtitle("")  title("Employment rate, by age", size(medsmall)) legend(order(2 "Age 22-30" 3 "Age 31-40") col(2)) ///
> saving("`output'/`dofile'_agegp2_empl.gph", replace)
file /accounts/projects/jr_ra/GRscarring/erratum/results/fig_ur_age_agegp2_empl.gph saved

. 
. gen altrecessionD=18*recession+74

. twoway area altrecessionD yearmo if agegp==3, yaxis(1) color(gs13) base(74) || ///
>  line empl_sa_m9 yearmo if agegp==3, yaxis(1) lstyle(p1) || ///
>  line empl_sa_m9 yearmo if agegp==4, yaxis(1) lstyle(p2) lpattern(longdash) || ///
>  line empl_sa_m9 yearmo if agegp==2, yaxis(1) lstyle(p3) lpattern(shortdash) || ///
>  , ylabel(74 (2) 92) ytitle("Employment rate") ///
>    xlabel(246 "1980" 306 "1985" 366 "1990" 426 "1995" 486 "2000" 546 "2005" 606 "2010" 666 "2015" 726 "2020") ///
>    xtitle("")  title("Employment rate of young college graduates, by age" , size(medsmall)) legend(order(2 " Age 22-25 " 3 " Age 2
> 6-29" 4 " Age 30-40") col(3)) ///
>    saving("`output'/`dofile'_agegp3_empl.gph", replace)
file /accounts/projects/jr_ra/GRscarring/erratum/results/fig_ur_age_agegp3_empl.gph saved

. 
. save `scratch'/`dofile'.dta, replace
(file /accounts/projects/jr_ra/GRscarring/erratum/scratch/fig_ur_age.dta not found)
file /accounts/projects/jr_ra/GRscarring/erratum/scratch/fig_ur_age.dta saved

. 
. 
. 
. ! rm `scratch'/extractcps.dta


. 
. 
. if `doasproject'==1 {
.       project, creates(`output'/`dofile'_base.gph)
project GRscar_erratum > do-file creates: "results/fig_ur_age_base.gph" filesig(1941232199:22113)
.       project, creates(`output'/`dofile'_agegp2_noempl.gph)
project GRscar_erratum > do-file creates: "results/fig_ur_age_agegp2_noempl.gph" filesig(690595442:23154)
.       project, creates(`output'/`dofile'_agegp2_unem.gph)
project GRscar_erratum > do-file creates: "results/fig_ur_age_agegp2_unem.gph" filesig(3256286019:23079)
.       project, creates(`output'/`dofile'_agegp2_empl.gph)
project GRscar_erratum > do-file creates: "results/fig_ur_age_agegp2_empl.gph" filesig(736807788:23116)
.       project, creates(`output'/`dofile'_agegp3_empl.gph)
project GRscar_erratum > do-file creates: "results/fig_ur_age_agegp3_empl.gph" filesig(3118633786:34659)
.        project, creates(`scratch'/`dofile'.dta)
project GRscar_erratum > do-file creates: "/scratch/public/jr_ra/GRscarring2024/erratum/scratch/fig_ur_age.dta" filesig(3306211396:4
> 47644)
.     }

. 
. 
end of do-file
      name:  plog_847
       log:  /accounts/projects/jr_ra/GRscarring/erratum/programs/analysis/fig_ur_age.log
  log type:  text
 closed on:  27 Nov 2024, 17:58:31
------------------------------------------------------------------------------------------------------------------------------------
